Kondo resonances and anomalous gate dependence of electronic conduction in 

single-molecule transistors 
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We report Kondo resonances in the conduction of single-molecule transistors based on transition 
metal coordination complexes. We find Kondo temperatures in excess of 50 K, comparable to those 
in purely metallic systems. The observed gate dependence of the Kondo temperature is inconsistent 
with observations in semiconductor quantum dots and a simple single-dot-level model. We discuss 
possible explanations of this effect, in light of electronic structure calculations. 



PACS numbers: 73.22.-f,73.23.Hk,85.65.+h 

In the Kondo HamiltonianQ, one of the most well- 
studied many-body problems in physics, an unpaired spin 
localized in a singly occupied electronic level is coupled 
via tunneling to an electronic bath. On-site Coulomb re- 
pulsion forbids real double occupancy of the level, but 
virtual processes favor antiferromagnetic exchange be- 
tween the local spin and the electronic bath. As T 
is reduced below a characteristic Kondo temperature, 
Tki these exchange processes "screen" the local mo- 
ment. The Kondo problem has undergone a resurgence, 
with atomic-scale studies of Kondo physics by scanning 
tunneling microscopy (STM)0, H and the realization 
of tunable Kondo systems in semiconductor quantum 
dots0, H, IE 0- With the recent development of single- 
molecule transistors (SMTs) based on individual small 
molecules 0, Kondo systems now include organometal- 
lic compound s El Ho) and fullerenes with normal^]] and 
ferromagnetic |12{ leads. 

In this Letter, we report Kondo physics in SMTs incor- 
porating transition metal complexes designed to contain 
unpaired electrons. As a function of gate voltage, Vq, we 
observe transitions from Coulomb blockade conduction 
to Kondo conduction, manifested as a strong peak in the 
differential conductance, G = dln/OVsr,, at zero bias, 
Vsd = in one charge state. At fixed Vq, the temper- 
ature dependence of the conductance peaks' amplitudes 
and widths agree well with the expected forms for spin- 
1/2 Kondo resonances. Observed SMT Kondo tempera- 
tures are >^ 50 K, comparable to those in purely metallic 
Kondo systems. We find that Tk(Vg) is strongly incon- 
sistent with the simple model of Kondo physics as seen in 
semiconductor devices 0,0- We discuss explanations for 
this anomalous gate dependence in light of spin-resolved 
electronic structure calculations of the complexes. 

Figure shows the simple single-level picture often 
used to describe Kondo physics in single-electron tran- 
sistors (SETs) . The intrinsic width of the single particle 
level is r = Ts + Td, determined by overlap of the single- 



particle state with the conduction electron states of the 
source and drain. The charging energy E c is the Coulomb 
cost of adding an extra electron to the molecule. The 
energy difference between the singly occupied level and 
the source/drain chemical potential is e, which is zero at 
charge degeneracy and varies linearly with Vq- 



(a) 





u 




Drain 



(b) NCS .SCN 





-30 30 
60 V SD [mV] 

V [V] (d) 



FIG. 1: (a) Energy level diagram of a typical SET in the 
Kondo regime, (b) Structural formula of transition metal 
complex before self-assembly, (c) Map of G(Vsr>, Vg) showing 
transition to Kondo resonance at 5 K. Brightness scales from 
G = (black) to 0.3 x 2e 2 /h. (d) G vs. Vsd at V G = 50 V 
(as indicated with arrow in (c)) at (top to bottom) 5 K, 16 K, 
and 30 K. 

Figure^ shows the structure of the neutral transition 
metal complexes measured in this study, as-synthesized. 
Two planar, conjugated ligands provide an octahedral co- 
ordination (compressed along the interligand (z) axis) of 
a transition metal ion, M. We have examined complexes 
with M = Co, Cu, and Zn, as well as individual ligands 
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and alkane chains. The Kondo devices in this Letter con- 
tain Co(II) and Cu(II). The as-synthesized molecules are 
characterized by EPR, SQUID, x-ray diffraction, cyclic 
voltammctry (CV), IR, and Raman spectroscopy. The 
Co(II) complex is high spin (3/2), with CV in solution 
confirming easily accessible Co(II)^Co(III) redox tran- 
sitions. The Cu(lf) complex is spin 1/2, with CV sup- 
porting easily accessible Cu(II)^->Cu(I) redox transitions. 
These complexes self-assemble on Au in tetrahydrofuran 
(THF) through loss of the -CN moieties and formation of 
Au-S covalent bonds0,0,^|. In principle, the ligands 
and remaining -SCN groups can also change their charge 
state. 

Devices are fabricated using Au source and drain 
electrodes on degenerately doped p+ silicon substrates 
which are used as the gate, with 200 nm of gate ox- 
ide. Source/drain electrodes are created by the controlled 
electromigration[Tfi| from lithographically defined Ti/Au 
(1 nm/15 nm thick) constrictions that are exposed for 
1 minute to oxygen plasma after liftoff, followed by self- 
assembly of molecules in solution (2 mM in THF) for 48 
hours. After rinsing and drying with dry nitrogen, sub- 
strates are placed in a variable temperature probe station 
for electromigration and measureme nt, with fabrication 
statistics comparable to past results flUllij . 

DC measurements of Id — Vsd are performed with the 
source electrode grounded, at various Vq. Differential 
conductance is computed by numerical differentiation, 
with spot checks by lock-in amplifier techniques. Device 
stability limits | Vsd | < 100 mV, while gate oxide limits 
| Vq I < 100 V. We consider only devices with significant 
gate response such that a charge transition is detected in 
maps of G(Vsd, Vg), to distinguish molecule-based effects 
from artifacts (e.g. metal nanoparticles). 

Figure^ shows a conductance map for a typical de- 
vice exhibiting a Kondo resonance. We have sufficient 
gate coupling to observe only a single charge degeneracy 
point for each device. From the slopes of the boundaries 
of the blockaded regime, one can estimate the constant 
of proportionality between changes in eVc and e as ex- 
pected from the model of Fig. QJi- For the device shown, 
Se » (CG/Ctot)eSVG » 10 -4 eSVo. The average value of 
this coefficient is 10~ 3 . The width of the charge degen- 
eracy resonance on the blockaded side of the degeneracy 
point sets an upper limit on T for the level participating 
in the redox state change. The T — ► limit of that width 
is proportional to T, while the degeneracy resonance is 
thermally broadened at finite temperature. Typical T 
values inferred from 5 K data in these devices are 3- 
30 meV. We note that the edges of Coulomb blockade 
diamonds, while usually distinct in the non-Kondo charge 
state, are often much weaker or apparently absent in the 
Kondo charge state, as also mentioned in Ref. [l2| . 

Fig. shows G(Vsd,Vg = 50 V) at three tempera- 
tures. The resonant peak decreases in magnitude while 
increasing in width as T is increased. Figs. |2K,b show 
the peak height and width, respectively, as a function 
of temperature. The solid line in (a) is a fit to the 



semiempirical expression |4| for the spin- 1/2 Kondo reso- 
nance in conduction, G(T) = GM + 2 1 / s ~ 1 ^V ;i , with 
s = 0.22. After the subtraction of a smooth background 
conductance, the adjustable parameters are the overall 
conductance scale, G c , and Tk. Similarly, the solid line 
in (b) is a fit to the expected form[3( for the full-width 
at half-maximum, FWHM = ^ {nk B T) 2 + 2(k B T K ) 2 , 
where the only adjustable parameter is Tk. The numer- 
ical coefficients are accurate to within a factor of order 
unity; other observations 0, E3 find better quantitative 
consistency with Tk as inferred from G(Vsd = 0, T) if 
the FWHM is assumed to be w (2fc B T K )/e. We adopt 
this latter assumption, and find good quantitative consis- 
tency between Tk values inferred from resonance height 
and width. 
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FIG. 2: (a) Temperature dependence of Kondo resonant peak 
height for the device of Fig.0; at Vg = 50 V. Solid line is the 
expected semiempirical functional form for spin-1/2 Kondo, 
with Tk = 69 K. (b) Temperature dependence of Kondo peak 
FWHM for the same device and Vg, with fit to expected 
functional form for spin-1/2 Kondo. Setting FWHM in the 
low-T limit to 2fcBTk/e gives Tk = 65 K. (c,d) Histograms 
of Tk as inferred from peak widths for Co and Cu complex 
devices, respectively. 

Figures^, d are histograms of Kondo temperatures in- 
ferred from resonance widths for 26 Co-containing Kondo 
devices (out of 921 electrode pairs examined at low tem- 
peratures), and 12 Cu-containing Kondo devices (out of 
397 electrode pairs examined at low temperatures) . Only 
devices exhibiting zero-bias resonances and clear charge 
degeneracy points are considered here. The observed 
weak dependence of Tk on Vq (see Fig. |3] and later dis- 
cussion) means that these distributions are relatively in- 
sensitive to the choice of Vg at which Tk is inferred. No 
Kondo resonances were observed in 370 control devices 
using alkanethiol chains, bare metal electrodes, and elec- 
trodes exposed to solvents and poor vacuum. 

Kondo physics in these complexes is clearly strong, 
with Tk values similar to those reported in STM mea- 
surements of Co atoms on Au(lll)|3- The Kondo tem- 



perature is expected to be 0,0: 
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outside the mixed valence (e/T < 1) regime. For small 
molecules, E c is likely to be hundreds of meV, while T is 
empirically tens of meV. The prefactor in Q then implies 
Tk can be as large as hundreds of Kelvin. 




FIG. 3: Log T K (normalized) inferred from FWHM of 5 K 
Kondo peak in G(Vsd), as a function of Vg normalized by 
the width of the Coulomb blockade charge degeneracy point 
(= e/T in the simple model of Fig.^) for several devices. In 
the model of Fig. [TJu, a parabolic dependence of log Tk on Vg 
is expected. Values of T inferred for these devices are, top 
down, 22.6, 11.5, 18, 3.3, 12.6, 14.2, and 26.8 meV. 

Assuming a model as in Fig. e = at the charge 
degeneracy point. Normalizing the difference in Vg away 
from charge degeneracy by the width in Vq of the charge 
degeneracy resonance in the Coulomb blockade regime 
gives a lower limit on e/T. Figure |21 shows (normalized) 
Tk as inferred from low-T resonance FWHM ~ 2fcB?k/e 
as a function of inferred e/T for several devices. Eq. (JTJ 
predicts a quadratic dependence of logTk on e, with a 
minimum in Tk at e = E c /2. The measured dependence 
of inferred Tk(Vq) ^ s much less steep; indeed, for Sample 
CuM7, Tk actually increases as Vq is shifted away from 
the nearest charge degeneracy. Thus, the simple model 
of Fig. which works well for semiconductor quantum 
dot experiments 0, 0, is inconsistent with this data. We 
note that the Tk(Vg) dependence reported in Ref. 9] for 
e/T > 1 is also surprisingly weak. 

We consider explanations for this deviation from sim- 
ple expectations. We dismiss as unlikely that the zero- 
bias peak may not be a true Kondo resonance, given (a) 
quantitative consistency of the functional forms for the 
resonance G(Vsu,T) with Kondo expectations; (b) the 
appearance of the zero-bias resonance coincident with 
passing through charge degeneracy points; (c) the lack 
of such resonances in control devices; (d) the similarity 
to other Kondo data reported in SMTs|9lllfJ|. 

The natural explanation for the anomalous gate depen- 
dence is that the normalized SVg used as e/T in Fig. is 
not the true e /T relevant to the Kondo Hamiltonian that 
gives Eq. In semiconductor quantum dots 6] infer- 
ring e/T by normalizing SVq is quantitatively consistent 



with Eq. Presumably some mechanism intrinsic to 
the molecular system renormalizes either the effective e, 
the effective T, or both, away from the simple picture of 
Fig.DJi.. 

Orbital degeneracy for the unpaired spin is one 
possibility^. When an TV-fold degeneracy exists, Tk 
as defined by Eq. is enhanced, such that the denomi- 
nator of the exponent becomes NTE C . Thus the normal- 
ization of the abcissa in Fig. [3] would effectively be too 
large by a factor of N, qualitatively explaining the appar- 
ently weak dependence of Tk{Vq) if N ~ 5. However, 
EPR spectra of both the Co and Cu-based complexes 
in solution phase show no indication of such a large de- 
generacy. Furthermore, we have performed spin-resolved 
calculations on the transition metal complexes to exam- 
ine their electronic structure, as shown in Fig. 01 (further 
details are available |lj|). In neither complex are large 
degeneracies expected. While it is conceivable that the 
self- assembled compounds could have different orbital de- 
generacies than the isolated molecules, it seems unlikely 
that both complexes, with their differing isolated elec- 
tronic structures, would have such similar properties. 

The electronic structure calculations reveal an addi- 
tional energy scale that is often small in semiconductor 
dots, but in the molecular system is comparable ener- 
getically to the inferred T, E c , e, and expected single- 
particle level spacing: intramolecular exchange. In both 
Co and Cu complexes, intramolecular exchange is strong. 
In the Co case, to a good approximation ligand field ef- 
fects split the d-states into one doubly degenerate, one 
nearly doubly degenerate, and one singly degenerate set 
of states. The complex has a total spin of S = 3/2, as 
confirmed by EPR. A minority spin half-occupied d state 
is at the Fermi level with two nearly degenerate d states 
1.1 eV above the Fermi level and the remaining occupied 
minority d electron 0.5 eV below the Fermi level. Ex- 
change splittings pull the majority d-band down signifi- 
cantly, with the highest occupied majority d state 1.8 eV 
below the Fermi level. The fermiology is complicated by 
an additional delocalized unpolarized two-fold degener- 
ate carbon 2p state lying 0.05 eV below the Fermi level. 
Also there is a delocalized doubly degenerate molecular 
state 0.85 eV below the Fermi level. The Cu complex has 
a total spin of S — 1/2, with the Fermi level unchanged 
in comparison to Co. This means that the minority spin 
d states are pulled to lower energies by approximately 
1.1 eV as shown in the figure. The structure calcula- 
tions also confirm significant derealization of the major- 
ity spin states, extending into the ligands, consistent with 
large T values for these systems. Jahn- Teller distortion 
is not considered, since the substrate-complex interaction 
is expected to be large compared to this effect. 

The occasional observation of Kondo resonances with 
similar Tk values on both sides of charge degeneracy 
points further indicates that intramolecular exchange 
cannot be neglected in these devices. However, it is not 
immediately clear how this strong interaction could lead 
to the observed phenomenon. 
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FIG. 4: Spin-resolved projected densities of states of the 
Cu(II) and Co(II) complexes, respectively. Minority (upper) 
and majority (lower) spin levels are indicated, demonstrat- 
ing the strong effects of intramolecular exchange. The black 
curve is the total DOS, and the dashed (blue) is the DOS 
projected onto the transition metal d states. Further details 
are available p^l . 

Vibrational effects are another piece of physics in- 
trinsic to the molecular system. The signature of 
electron-phonon interactions has been observed in inelas- 
tic electron tunneling spectra of SMTs made with these 
complexes 113 . M oderate coupling of charge and vibra- 
tional modesjl^ localized to the molecule can strongly 
both increase Tk and decrease its gate dependence rela- 
tive to the case with no vibrational coupling |20|] . A quan- 
titative estimate of electron- vibrational couplings would 
facilitate testing this hypothesis, and should be obtain- 
able from further quantum chemistry calculations. 

Finally, it is also possible [2l| that screening correla- 



tions in the mixed valence regime can renormalize the 
measured T to a value different than the T relevant to 
the Kondo temperature. Full quantum chemistry calcu- 
lations of molecules bound to realistic Au leads includ- 
ing Kondo and many-body correlations (as done for Co 
atoms on Au(lll)[22|) are essential to a better under- 
standing the observed effects, and are beyond the scope 
of this paper. 

In measuring the electronic properties of single- 
molecule transistors containing transition metal com- 
plexes, we observe strong Kondo physics, indicating that 
conjugated ligands can provide extremely effective cou- 
pling of spin degrees of freedom to metal leads. We also 
find a Tk(Vq) that is vastly weaker than that seen in 
semiconductor quantum dot realizations of the Kondo 
effect and expected for the simple model of Fig. We 
have discussed possible explanations for this anomalous 
dependence in light of electronic structure calculations 
of the complexes. While a complete understanding will 
require more sophisticated modeling and further mea- 
surements, these data demonstrate that correlated states 
involving SMTs can exhibit rich effects not seen in their 
semiconductor quantum dot counterparts. 
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